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' - ' Abstract. In this article we show how to compute a matrix representation and 

^—^ the implicit equation by means of the method developed in IBotlObI , using the 

^^ computer algebra system Macaulay2 |GS|. As it is probably the most interesting 

1—^ case from a practical point of view, we restrict our computations to parametriza- 

, ^ tions of bigraded surfaces. This implementation allows to compute small exam- 

' pies for the better understanding of the theory developed in IBotlObI , and is a 

T— H complement to the algorithm IBDIOI . 
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1. Introduction and background 



The interest in computing explicit formulas for resultants and discriminants 
goes back to Bezout, Cayley, Sylvester and many others in the eighteenth and 
Ch nineteenth centuries. The last few decades have yielded a rise of interest in the 

implicitization problem for geometric objects motivated by applications in com- 
^-H puter aided geometric design and geometric modeling as can be seen in for exam- 

> pie in ||Kal911 IMC92b[ IMC92ai IAGR95I ISC95l . This phenomena has been empha- 

O sized in the latest years due to the increase of computing power, see for example 

^ llASOTl ICo^^ IBCD031 [BJ031 iBCOSl |BCJ09| IBD071 IBot09l IBDD091 lEBtTOal iBotlOb J . 

-y-s Under suitable h3rpotheses, resultants give the answer to many problems in 

• elimination theory, including the implicitization of rational maps. In turn, both 

L». resultants and discriminants can be seen as the implicit equation of a suitable 

(^ map (cf. IIDFS07I ). Lately, rational maps appeared in computer-engineering con- 

T-H texts, mostly applied to shape modeling using computer-aided design methods 

L/ for curves and surfaces. A very common approach is to write the implicit equa- 

• ^H tion as the determinant of a matrix whose entries are easy to compute. In gen- 

/\^ eral, the search of formulas for implicitization rational surfaces with base points 

is a very active area of research due to the fact that, in practical industrial de- 
sign, base points show up quite frequently. In IIMC92al , a perturbation is applied 
to resultants in order to obtain a nonzero multiple of the implicit equation. In 
l lBJOSl JBCOS. .BCJ091 iBDOTl IBDD09I iBotTOal IBotlObI show how to compute the im- 
plicit equation as the determinant of the approximation complexes. 

In HBotlObl we present a method for computing the implicit equation of a hyper- 
surface given as the image of a rational map (p : ^ ---* P", w^here ^ is an arith- 
metically Cohen-Macaulay toric variety. In IIBDD09II and HBotlOal , the approach 
consisted in embedding the space ^ in a projective space, via a toric embedding. 
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The need of the embedding comes from the necessity of a Z-grading in the coordi- 
nate ring of !X , in order to study its regularity. We exploit the natural structure of 
the homogeneous coordinate ring of the toric variety where the map is defined. 

In [BotlObJ we introduce the "good" region in G where the approximation com- 
plex Z, and the symmetric algebra Sym^(/) has no S-torsion. Indeed, we define 
for 7 e G, 

^b(7):= U (6s(7)-fc-7)cG, 

0<fe<min{m,cdi3(_R,)} 

where6s(7) := Ufc>o(SuPPG(^B(-f^)) + ^'7)- This goes in the direction of proving 
the main theorem loc. cit., IIBotlObI Thm. 4.4 and Rmk. 4.5]. Precisely it asserts 
that, w^hen ^ is a (n — l)-dimensional non-degenerate toric variety over a field K, 
and S its Cox ring (cf. IICox95ll '). For a rational map : J^" ---> P" defined by n + 1 
homogeneous elements of degree p e C\{^). If dim(y(/)) < in JT and V{I) 
is almost a local complete intersection off V{B), we proved in Theorem [iBotlObI 
Thm. 4.4] that, 

det((Z.)^) = i/deg(0) . Q £ ]K[T], 

for all 7 ^ ^b{p), w^here H stands for the irreducible implicit equation of the 
image of (p, and G is relatively prime pol5momial in K[T]. This result is a new 
formulation of that in pj03| Thm. 5.7] and BBotlOai Thm. 10 and Cor. 12] in this 
new setting. 

In this article we show how to compute a matrix representation and the unplicit 
equation with the method developed in BBotlObl , following BBotlOal , using the 
computer algebra system Macaulay2 [GS]. As it is probably the most interesting 
case from a practical point of view, we restrict our computations to parametriza- 
tions of a bigraded surface given as the image of a rational map : P^ x P^ -— > P-^ 
given by 4 homogeneous polinomials of bidegree (£1,62) € Z^. Thus, in this 
case, the Z-complex can be easyly computed, and the region 9^^(61,62) where 
it is acyclic is 

^5(61,62) = (SuppG(il|(i?)) + (61,62)) U (SuppG(JT|(i?)) +2 • (61,62)). 

This implementation allows to compute small examples for the better under- 
standing of the theory, but we are not claiming that this implementation is opti- 
mized for efficiency; anyone trying to implement the method to solve computa- 
tionally involved examples is weU-advised to give more ample consideration to 
this issue. 

2. Implementation in Macaulay 2 
Consider the rational map 

/J^ Pl X Pl --^ P3 

(s : u) X (i : «) K^ (/i : /2 : /3 : h) 

where the polynomials fi = fi{s, u, t, v) are bihomogeneous of bidegree (2, 3) e Z^ 
given by 

• /i = s'^t^ + "^sut^ + iuH'^ + As^t^v + Ssut^w + Qu'^t^v + Is^tv"^ + Ssutv"^ + 

• /2 = 252^3 _ 3s2^2^ „ g2^^2 _(_ g^£2^ ^ Ssutv'^ ~ Su^t^w + 2u'^tv^ - U^v'^ , 



AN ALGORITHM FOR COMR IMPL. EQ. OF BIGRADED RAT. SURFACES 

fi = is^fv - 2sut^ - s^tv^ + sut^v - 3sutv^ - uH'^v + Au^tv"^ - v?v'^. 



Our aim is to get the implicit equation of the hypersurface im(/) of P'^. 
First we load the package "Maximal minors" 

il : load "maxminor .m2" 

Let us start by defining the parametrization / given by (/i, /2, /s, A)- 
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S=QQ[s,u,t, v,Degrees=>{ {1,1,0}, {1,1,0}, {1,0,1}, {1,0,1}}: 

el=2; 

e2 = 3; 



15 : fl = l*s"2*t"3 + 2*s*u*t"3 + 3*u~2*t'3 + 4*s"2*t''2*v+5*s*u*t"2*v+6*u"2*t"2*v+ 

7*s"2*t*v"2+8*s*u*t*v"2+9*u'2*t*v"2+10*s"2*v"3+l*s*u*v'3+2*u"2*v~3; 

16 : f2^2*s"2*t'"3-3*s"2*t"2*v-s'"2*t*v"2 + s*u*t'"2*v+3*s*u*t*v'"2-3*u"2*t'"2*v+ 

2*u"2*t*v'"2-u"2*v"3; 

17 : f3^2*s "2*t "3-3*3 "2*t"2*v-2*s*u*t" 3 + s '"2*t*v'"2 + 5*s*u*t '"2*v-3*s*u*t*v'"2- 

3*u"2*t'"2*v+4*u"2*t*v'"2-u'"2*v"3; 

18 : f4 = 3*s"2*t'"2*v-2*s*u*t"3-s'"2*t*v"2 + s*u*t'"2*v-3*s*u*t*v'"2-u"2*t"2*v+ 

4*u"2*t*v'"2-u"2*v"3; 

We construct the matrix associated to the pol5niomials and we relabel them in 
order to be able to automatize some procedures. 

19 : F=matrlx{ {fl, 12, f3, f4} } 

o9 = I S2t3+2sut3+3u2t3+4s2t2v+5sut2v+6u2t2v+7s2tv2+8sutv2+9u2tv2+10s2v3+ 

suv3+2u2v3 2s2t3-3s2t2v+sut2v-3u2t2v-s2tv2+3sutv2+2u2tv2-u2v3 

2s2t3-2sut3-3s2t2v+5sut2v-3u2t2v+s2tv2-3sutv2+4u2tv2-u2v3 

-2sut3+3s2t2v+sut2v-u2t2v-s2tv2-3sutv2+4u2tv2-u2v3 

1 4 
o9 : Matrix S < S 

The reader can experiment with the implementation simply by changing the 
definition of the pol3niomials and their degrees, the rest of the code being identical. 

Let A: be a field. Assume ,9^ is the biprojective space Pj^ x Pj|. Take Ri := 
k[xi,X2],R2 ■■= k[yi,y2,y3,y4],andG := 1? . Write i? := i?i«)fe-R2 andsetdcg(a:;i) = 
(1,0) and deg(2/j) = (0,1) for all i. Set ai := {xx,X2), 02 := (yi, ^2,2/3, 2/4) and 
define i? := ai • 02 C i? the irrelevant ideal of R, and m := ai + 02 C i?, the ideal 
corresponding to the origin in Spec(i?). 

By means of the Mayer- Vietoris long exact sequence in local cohomology, we 
have that: 

(1) lil{R)^u,\®^ujl^, 

(2) Hl{R)=Rl,{R)=^l, 

(3) B^s (R) = for alH 7^ 2 and 3. 

Thus, we get that: 

(1) Suppc(i/|(i?)) = -N X N + (-2, 0) U N X -N + (0, -2). 

(2) SnppciHUR)) = -N X -N + (-2, -2), . 
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(-2,0) 



• 

(-1,-1) 



(0,0) 




And thus, 

$Kb(2, 3) = {SuppciHliR)) + (2, 3)) U (SuppG(i/3 (i?)) + 2 • (2, 3)). 
Obtaining 

C<nB(2, 3) ^ (N^ + (1, 5)) U (N^ + (3, 2)). 
As we can see in Example 12} a Macaulay2 computation gives exactly this region 
(illustrated below) as the acyclicity region for Z, . 

ilO : nu={5,3,2}; 

An alternative consists in taking 

ilO : nu={6,l,5}; 




Non-vanishing of 
local cohomology 



Anyhow, it is interesting to test what happens in different bidegrees i/ e Z^ by 
just replacing the desired degree in the code. 
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In the following, we construct the matrix representation M. For simplicity, we 
compute the whole module Zi, which is not necessary as we only need the graded 
part (Zi)^Q. In complicated examples, one should compute only this graded part 
by directly solving a linear system in degree vq. 



ill 
il2 
113 
114 



ZO=S'l; 

Zl=kernel koszul(l,F); 
Z2^kernel koszul(2,F); 
Z3^kernel koszul(3,F); 



115 : d={el+e2,el,e2} 

116 : hfZOnu = hllbertPunctlon (nu, ZO ) 

016 = 12 

117 : hfZlnu = hllbertPunctlon (nu+d, Zl ) 

017 = 12 

118 : hfZ2nu = hllbertFunction (nu+2*d, Z2 ) 

018 = 

119 : hfZ3nu = hllbertFunction (nu+3*d, Z3 ) 

019 = 

120 : hfnu = hf ZOnu-hf Zlnu+hf Z2nu-hf Z3nu 

020 = 

Thus, when i^o = (3, 2) or vq = (1, 5), we get a complex 

M 

(Z.)^„ : ^ ^ ^ K[X]i2 _^ ]f^[x]i2 ^ 0. 

and, hence, dct((Z,)j^j,) = dct(M^Q) e IK[X]i2 is an homogeneous polynomial of 
degree 12 that vanishes on the closed image of 0. We compute here the degree of 
the MacRae's invariant which gives the degree of dci{{Z,)v^). 

121 ihilbertFunctlon (nu+d, Zl) -2*hilbertFunctlon (nu+2*d, Z2) + 

3*hilbertFunctlon (nu+3*d, Z3) 

021 = 12 

122 : GG=ideal F 

23 32322 2 22 22 

022 ^ Ideal (s t +2s*u*t +3u t +4s t v+5s*u*t v+5u t v+7s t*v + 

222 23 323 2322 2 

8s*u*t*v +9u t*v +10s V +s*u*v *2u v , 2s t -3s t v+s*u*t v- 

2222 22223 23 322 

3u t v-s t*v +3s*u*t*v +2u t*v -u v , 2s t -2s*u*t -3s t v+ 

22222 22223 3 

5s*u*t v-3u t v+s t*v -3s*u*t*v +4u t*v -u v , -2s*u*t + 

22 22222 22223 

3s t v+s*u*t v-u t v-s t*v -3s*u*t*v +4u t*v -u v ) 

o22 : Ideal of S 

123 : GGsat=saturate (GG, Ideal (s, t) *ideal (u, v) ) 
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22 222 22 22223 

o23 ^ ideal (3s t v-3s*u*t v-u t v-3s t*v +3s*u*t*v +2u t*v -u v , 

23 2 22 22 22223 

9u t +42s*u*t v+28u t v+45s t*v -15s*u*t*v +19u t*v +30s v + 

323 3 222 22223 

3s*u*v +13u V , s*u*t -2s*U'^t v-s t*v +3s*u*t*v -u t*v , s t - 

2 22 22 22223 424 

s*u*t v-2u t v-2s t*v +3s*u*t*v +2u t*v -u v , 30s*u*v -u v , 

24 242324 324 23 24 

15s V +14u V , u t*v -u V , 30s*u*t*v -u v , 15s t*v +14u v , 

2 2 2 2 4 

U t V -U V ) 

023 : Ideal of S 

124 : degrees gens GGsat 

024 = {{{0, 0, 0}}, {{5, 2, 3}, {5, 2, 3}, {5, 2, 3}, {5, 2, 3}, {6, 

2, 4}, {6,2, 4}, (5, 2, 4}, {6, 2, 4}, {6, 2, 4}, (6, 2, 4}}} 

024 : List 

125 : H=GGsat/GG 

025 ^ subquotlent (I 3s2t2v-3sut2v-u2t2v-3s2tv2+3sutv2+2u2tv2-u2v3 

9u2t3+4 2sut2v+2 8u2t2v+4 5s2tv2-15sutv2+19u2tv2+3 0s2v3+3suv3+ 
13u2v3 Sut3-2sut2v-s2tv2+3sutv2-u2tv2 S2t3-sut2v-2u2t2v-2s2tv2+ 
3sutv2+2u2tv2-u2v3 30suv4-u2v4 15s2v4+14u2v4 u2tv3-u2v4 
30sutv3-u2v4 15s2tv3+14u2v4 u2t2v2-u2v4 j, | s2t3+2sut3+3u2t3+ 
4s2t2v+5sut2v+6u2t2v+7s2tv2+8sutv2+9u2tv2+10s2v3+suv3+2u2v3 
2s2t3-3s2t2v+sut2v-3u2t2v-s2tv2+3sutv2+2u2tv2-u2v3 2s2t3-2sut3- 
3s2t2v+5sut2v-3u2t2v+s2tv2-3sutv2+4u2tv2-u2v3 -2sut3+3s2t2v+ 
Sut2v-u2t2v-s2tv2-3sutv2+4u2tv2-u2v3 | ) 

1 

025 : S~module, subquotlent of S 

126 : degrees gens H 

026 = {{{0, 0, 0}}, {{5, 2, 3}, {5, 2, 3}, (5, 2, 3}, {5, 2, 3}, {6, 

2, 4}, {6,2, 4}, {6, 2, 4}, {6, 2, 4}, {6, 2, 4}, {6, 2, 4}}} 
o26 : List 

Now, we focus on the computation of the implicit equation as the determinant 
of the right-most map. Precisely, we will build-up this map, and later extract a 
maximal minor for taking its determinant. It is clear that in general it is not the 
determinant of the approximation complex in degree v, but a multiple of it. We 
could get the correct equation by taking several maximal minors and considering 
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the gcd of the determinants. This procedure is much more expensive, hence, we 
avoid it. 

Thus, first, we compute the right-most map of the approximation complex in 
degree u 

ill : R=S [T1,T2, T3, T4] ; 

128 : G=sub(F,R); 

1 4 

028 : Matrix R < R 

We compute a matrix presentation for (Zi)i, in Ki: 

129 :Zlnu=super basis (nu+d, Zl) ; 

4 12 

029 : Matrix S < S 

130 : Tnu=matrix{ {T1,T2,T3,T4} }*substitute (Zlnu,R) ; 

1 12 

030 : Matrix R < R 

131 : lll^matrlx {{s,t,u,v}} 

031 = 1 s t u V 

1 4 

031 : Matrix S < S 

132 : lll=sub(lll,R) 

032 = 1 s t u V 

1 4 
o32 : Matrix R < R 

133 : 11={}; 

134 : for 1 from to 3 do { ll=append (11, 111_ ( 0, 1) ) } 

Now, we compute the matrix of the map {Zi)iy -^ A^[Ti,T2,T3, T4] 

135 : (m,M) =coef f Icients (Tnu, Variables=>ll,Monomials^>substitute ( 

basis (nu, S) , R) ) ; 

13 6 : M; 

12 12 

o36 : Matrix R < R 

137 : T=QQ[T1,T2, T3, T4] ; 

138 : ListofTandO = { Tl , T2, T3, T4 } ; 

139 : for 1 from to 3 do { LlstofTandO=append (ListofTandO , ) } ; 

140 : p=map (T,R, ListofTandO) 

o40 = map(T,R, {Tl, T2, T3, T4, 0, 0, 0, 0}) 
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o4 : RingMap T < R 

141 :N=MaxCol (p(M) ) ; 

12 12 

o41 : Matrix T < T 

The matrix M is the desired matrix representation of the surface y. We can 
continue by computing the implicit equation by taking determinant. As we men- 
tioned, this is fairly more costly. If we take determinant what we get is a multiple 
of the implicit equation. One wise way for recognizing which of them is the im- 
plicit equation is substituting a few points of the surface, and verifying which 
vanishes. 

Precisely, here there is a multiple of the implicit equation (by taking several 
minors we erase extra factors): 

142 :Eq=det(N); factor Eq; 

We verify the result by substituting on the computed equation, the pol5momials 
/i to /4. We verify that in this case, this is the implicit equation: 

144 : use R; Eq=sub (Eq, R) ; 

146 : sub(Eq, {T1=>G_(0, 0) ,T2 = >G_(0, 1) ,T3=>G_(0,2) ,T4=>G_(0, 3) }) 

o46 = 

o46 : R 
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